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Abstract: We use the holographic AdS/QCD soft-wall model to investigate the spectrum 
of scalar glueballs in a finite temperature plasma. In this model, glueballs are described 
by a massless scalar field in an AdSs black hole with a dilaton soft-wall background. Using 
AdS/CFT prescriptions, we compute the boundary retarded Green's function. The cor- 
responding thermal spectral function shows quasiparticle peaks at low temperatures. We 
also compute the quasinormal modes of the scalar field in the soft-wall black hole geometry. 
The temperature and momentum dependences of these modes are analyzed. The positions 
and widths of the peaks of the spectral function are related to the frequencies of the quasi- 
normal modes. Our numerical results are found employing the power series method and 
the computation of Breit-Wigner resonances. 
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1. Introduction 

In the last years the AdS/CFT correspondence [1] has been recognized as an important 
tool to investigate diverse aspects of strongly coupled gauge theories. This correspondence 
relates a ten-dimensional string theory or eleven-dimensional M-theory to a conformal 
gauge theory on the corresponding spacetime boundary. The most celebrated version of 
this correspondence relates string theory in AdSs x space to four-dimensional SU(iVc) 
Yang- Mills gauge theory with large Nc and extended M = A supersymmetry. 

At low energy string theory is represented by an effective supergravity theory. Hence, 
the AdS / CFT correspondence implies a gauge / gravity duality. Using this duality, one can 
calculate quantities like quantum correlation functions for gauge-field boundary operators 
from the classical supergravity action for bulk fields [2,3]. 
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One can formulate a finite temperature version of the AdS/CFT correspondence in- 
cluding a black hole in the AdS spacetime. The Hawking temperature of the black hole, 
related to the horizon position, is identified with the temperature of the dual gauge the- 
ory [4]. Following this approach, one can calculate various thermal Green's functions for 
the boundary gauge theory by taking derivatives of the on-shell bulk action for fields sub- 
jected to certain conditions. In the Euclidean version of the AdS/CFT correspondence, 
the prescription to compute Green's functions is the same found in [2, 3] with the classi- 
cal bulk fields uniquely determined by their value at the boundary and the requirement 
of regularity at the horizon. In Minkowski space, a prescription for calculating retarded 
propagators at finite temperature was proposed in [5] (see also Refs. [6-8]) imposing that 
the fields satisfy purely incoming-wave condition at the horizon. Physically, this means 
that the black holes only absorb without any emission. An important consequence of this 
prescription is that the poles of the retarded correlation function of an operator in the 
boundary conformal field theory correspond to the quasinormal (QN) frequencies of the 
holographic dual field in the AdS black-hole spacetime. In the context of the AdS/GFT 
correspondence, quasinormal modes have been object of study of many works (see Ref. [9] 
for a recent review) . 

The AdS/CFT correspondence can be also an important tool to study many problems 
in strong interactions, like the calculation of hadronic spectra. The theory that describes 
strong interactions is QCD, a SU(3) Yang-Mills theory with a coupling constant that runs 
with the energy. At low energies the coupling is strong and one can not use perturbative 
methods. In this regime, one needs other approaches like lattice QCD or phenomenological 
effective models. Recently, the idea of modifying the AdS space introducing an infrared 
cut off led to phenomenological models to describe strong interactions. This modification 
implies a mass gap in the boundary gauge theory that breaks the conformal symmetry 
present in the AdS/CFT correspondence. This kind of approach is known as AdS/QGD. 

The simplest AdS/QCD model, known as hard wall, introduces a mass gap via a cut 
off on the bulk geometry. In this way it is possible to use the gauge/gravity duality to 
estimate hadronic masses [10-15]. The hard wall is confining at low temperatures [16] 
and deconfining at high temperatures [17]. Another simple model of AdS/QCD, known 
as soft-wall, introduces a mass gap through a non-uniform background scalar field in the 
AdS spacetime. This soft-wall model presents linear Regge trajectories for vector mesons 
and glueballs [18-20]. Confinement/deconfinement thermal phase transition in these two 
models was studied in [21-24]. It was found in [22] that this corresponds to a gravitational 
Hawking-Page phase transition [25] between spaces with and without a black hole. This 
is a first order phase transition, where two different phases can coexist. The black-hole 
spacetime is thermally favored at high temperatures, while at temperatures below some 
critical model-dependent value, the thermally favored space is a pure AdS. 

In this article we study the spectrum of scalar glueballs at finite temperature using the 
soft- wall AdS/QCD model. This model, at high temperatures, consists of an AdS space 
with a nonuniform background dilaton field and a black hole. At intermediate and low 
temperatures, we consider, respectively, the (supercooled) metastable and unstable phases 
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of the plasma where the black hole is still present . 

Our motivation for studying the supercooled phase of the plasma is to gain some insight 
into the mechanism of glueball formation when the plasma cools down. We know that at 
the high temperatures reached when the plasma is formed, there are no stable hadronic 
states. On the other hand, at zero temperature, when the plasma disappears, there should 
be only colorless hadronic states with a well defined mass spectrum. These two separate 
regimes have been extensively studied using gauge/gravity duality. We want to investigate 
in this article how the transition between the high temperature deconfined phase and the 
zero temperature confined phase takes place. We think that this kind of analysis may be 
useful to understand hadronic formation in heavy ion collisions. 

It is important to remark that the five dimensional background of the soft-wall phe- 
nomenological model does not arise as a solution of the 5D Einstein equations. In this 
model there is a time-independent background scalar field <I> that is a function of an in- 
frared energy scale. In this article we will study the dynamics of an extra scalar field 
in this background at finite temperature. As the infrared energy scale goes to zero, the 
background $ of the soft- wall model vanishes. In this limit, we recover the usual AdS/CFT 
picture and the scalar field (j) can be thought as a perturbation of the trivial dilaton vacuum 
or of the scalar part of the metric. In the AdS/CFT correspondence this field is dual to 
the scalar glueball operator. Motivated by this fact, we assume in this work that the field 
(p is dual to scalar glueballs even in the presence of the infrared scale. 

Therefore, we identify the quasinormal modes (QNM) of the massless scalar field in 
the bulk as the spectrum of scalar glueballs in the boundary. We study the zero- and finite- 
temperature retarded Green's functions of scalar glueball operators both with numerical 
and analytical techniques. The spectral function is identified with the imaginary part of 
the retarded Green's function. We perform also a numerical analysis to obtain the real 
and imaginary parts of the quasinormal frequencies as functions of the temperature and 
the spatial momentum. In such analysis a good convergence is achieved using the power 
series method [27] at high temperatures and the Breit-Wigner ressonance method [28] at 
low temperatures. 

2. Glueballs in the soft- wall model at zero and finite temperature 

In this section we obtain the equation of motion of a scalar field in the soft-wall model 
at zero and finite temperature and transform this equation into a Schrodinger-like form. 
We also describe the wave functions that are relevant for the computation of the glueball 
retarded Green's function and analyze the corresponding effective potential. 

2.1 The soft- wall model at zero temperature 

According to the AdS/CFT correspondence, the particle states on the boundary field theory 
are dual to normalizable fields in the bulk spacetime. As a consequence of the conformal 
symmetry, the mass spectrum of these states is continous. In the AdS/QCD models the 

'^A similar situation was considered recently for the case of mesons in the D3-D7 model [26]. 
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conformal invariance is broken by the presence of an infrared cut-off turning the particle 
mass spectrum discrete. The soft-wall model of Ref. [18] was originally formulated at zero 
temperature in order to reproduce the approximate linear Regge trajectories for vector 
mesons. It consists of a five-dimensional AdS spacetime and a background dilaton-like 
field ^{z). The AdSs metric in Poincare coordinates is given by 



-dt' 



+ Y^^dx'f + dz"' 



(2.1) 



where A{z) = — ln(z/L) and L denotes the AdS curvature radius. The interaction of the 
background field ^{z) with the bulk fields is defined through the replacement of the action 
integrals: 



(2.2) 



The factor exp[— <I>(2;)] with ^{z) = cz^ introduces an effective infrared cut-off in the AdS 
space, with -^/c playing the role of a mass scale. This modifies the normalizability condition 
for bosonic fields in the bulk. The solutions with plane wave dependence exp(iA;-2;) become 
normalizable only for particular discrete values of /c^ = — m-^, with nonvanishing m„ and 
n = 0, 1, 2, ... . This implies a mass gap in the dual gauge theory. For instance, the spectrum 
of vector mesons is [18] 



4c(n + 1) 



(2.3) 



where n is the radial quantum number. 

Scalar glueball states at zero temperature in the soft-wall model were studied in [19] by 
considering a massless scalar field living in the AdS spacetime ( |2.1| ) with the same dilaton 
background field ^ijc) = cz^ . The corresponding glueball spectrum is 



4c(n + 2) . 



(2.4) 



The parameter ^fc , in this case, sets the scale of the glueball masses. Other recent results 
on glueballs in holographic models can be found in [29]. For a recent review on glueballs 
see [30]. 



2.2 The finite temperature case 

The soft- wall model at finite temperature is composed by the dilaton field ^ 
AdS black-hole spacetime with translationally invariant horizon: 

3 



ds' 



-f{z)df + Y^{dx'f + f{z) 



-'dz^ 



i=l 



cz^ and an 



(2.5) 



where /(z) = 1 — {zjzh)^. The coordinate z is defined in the range < z < z/i and z = 
corresponds to the boundary of the space. The parameter Zh indicates the position of the 
event horizon, which is related to the black-hole Hawking temperature T by z/^ = l/vrT. 
According to the holographic dictionary, T also represents the temperature of the boundary 
field theory. 
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Thermo dynamics 

Considering the Euclidean version of the metrics given by (2.1) and ( p.5| ) with compactified 
times, Herzog [22] has shown that, in the presence of the dilaton field $ = cz^, there is a 
thermal competition between these spaces, which is characterized by the ratio T = -kT j yfc. 
The regularized gravity action densities for the thermal AdS and AdS black hole spaces 
are given by [24] 

'3 



iV2 
^ 'c 

4^2 



7 /3 



+ 



iV2 

47r 



1 



+ 



1 



(2.6) 
(2.7) 



c^z,i?^(-cz^)+e— W - 

where Ei{u) = — J^^dtt^^e"^ and = /3/vr. At low temperatures the thermal AdS 
space is thermodynamically preferred compared to the AdS black-hole space, since it has 
a smaller gravity action. Then a first-order Hawking-Page type phase transition between 
these two spaces occurs at ~ 2.38644. This transition is interpreted in the boundary 
field theory as a confinement /deconfinement phase transition as can be seen by calculating 
the corresponding entropy change. At T > Tc the AdS black-hole space has a smaller 
gravity action so it is the dominant configuration. This space is dual to the plasma phase 
of a strongly coupled Yang-Mills theory with large number of colors. In the left panel of 
Fig. H, we show the difference of the regularized gravity actions densities, Sbh — SAds, that 
determines the confinement/deconfinement phase transition. 
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Figure 1: Left: Difference of regularized gravity action densities in the soft- wall model as a 
function of [ttTY /c. Right: Specific heat density for the AdS black hole with the soft-wall dilaton 
as a function of (7rT)^/c. The quantities Sbh — S'^ds and Cbh arc in units of N'^c^/'^ /A-k. 



In this article, we consider the black-hole spacetime ( |2.5| ) for all temperatures. For 
temperatures lower than the black hole is not stable and can be interpreted as a super- 
cooled Yang-Mills plasma. If the sign of the specific heat is positive (negative) it will be in 
a metastable (unstable) phase. From the regularized action density for the soft-wall AdS 
black hole (|2.7D , we obtain the specific heat density 



Cbh — P 



IT 



2 I 



4c 



+ 6 



(2.8) 
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In the limit c — > we recover the usual (positive) result for the black hole without dilaton. 
In our case, the presence of the dilaton implies that the specific heat is negative for < 
0.75, as shown in the right panel of Fig. ^. 

Equation of motion 

The dual of the scalar glueball operator Tr(i^^) is a massless scalar field (/) in the bulk. In 
the finite temperature soft-wall model this field is described by the action 



d'xV^e-''g''''d,,(t>dN(^, (2.9) 



4^10 



where <&(-z) = cz^ is again the soft- wall dilaton field and Qmn is the black-hole metric given 
by ( |2.5| ). The constants L and kiq are the anti-de Sitter radius and the ten-dimensional 
gravitational constant. The indices M, N run over 0, 1, 4 where (^u = 0, 3) are the 
Ad boundary coordinates and = z is the extra radial coordinate. 

As usual in the soft-wall model, we do not consider here the backreaction of the dilaton 
<I> on the metric. From the action (|2.9|), we find the scalar field (p satisfies the equation 



dz {V^e-^g^'d^cp) +g^'^d^d,(i> = 0. (2.10) 



Introducing the Fourier transform with respect to the coordinates x^, 

4>{z,x) = J ^e''-^4>i^,k), (2.11) 

equation ( 2.10| ) takes the form 

e^/(z)9, (e-^/(z)9,0) + (a;^ - f{z)q^) <^ = , (2.12) 
where k^, = {-uj,qi), = Y.^^^ qf and B = ^ - 3A = cz'^ + 31n(2;/L). 
2.3 The Schrodinger equation and the wave functions 

In the study of classical- field perturbations around black holes, it is usually convenient 
to introduce the Regge- Wheeler tortoise coordinate r*. This coordinate is defined by the 
relation dr^ = —f{z)dz with z in the interval < z < z^. Since in our case l/f{z) is 
integrable we find the expression 



1 



, Z \ 1 / Zh 

arctan — H — in 



ZhJ 2 \Zh + z 



(2.13) 



Note that the horizon z = z^ corresponds to —oo and the boundary z = is chosen at 

r^, = 0. It is also convenient to use a Bogoliubov transformation of the form ip = e~^/^(^. 
Then, Eq. ( 2.12| ) reduces to a one-dimensional Schrodinger equation: 

d'^^i; + LO^ip = ViJj, (2.14) 
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where the effective potential V is given by 



(2.15) 



This potential can be exphcitly written in terms of the coordinate z as 



V{z) 



zV + ^ + ^^ + 2c.^ (l + ^)+cV/(.] 



(2.16) 



By investigating the asymptotic behavior of the Schrodinger equation ( |2.14| ) near the 
boundary z = we find two particular real solutions Tpi and ip2- These solutions have the 
following asymptotic form: 



^"2 



■01 = 
-3/2 



,5/2 



[l + aiiz^ + ai2z'^ H ] , 

[l + 021^^ + 022^;'' H ] + & V"! ln(cz^ 



(2.17) 
(2.18) 



where the ellipses denote higher powers of z and the coefficients above are given by 



UJ 



ail 



0.21 



q'-2c 



12 



2c 



0-12 



b 



(w^ -q'^ - 2cf 



12 X 32 



+ 



1 
2^ 



4 + 32' 



4c) 



(2.19) 
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The wave function ij^i is normalizahle due to its asymptotic behavior z^/^, while the wave 
function ^/'2 is non-normalizahle since it behaves asymptotically as z~^/'^. The coefficient 
022 is arbitrary and, in particular, we can choose 022 = 0. The solutions ijji and ^2 will be 
useful to construct the glueball retarded Green's function in section ^ and to compute the 
black-hole quasinormal modes in section ^ 

Another special kind of solutions to the Schrodinger equation (2.14) are the ingoing 
and outgoing wave functions, ip- and '(/'+• These solutions are motivated by the fact that 
the potential V{z) vanishes when z — > z/^ so that Eq. ( |2.14| ) takes the form of a free 
particle equation with solutions proportional to exp(ibit(;r^,). Since in our convention for 
the Fourier transform the temporal dependence of (j) is exp(— iwf), the minus (plus) sign 
can be interpreted as an ingoing (outgoing) plane wave in the sense that the wave travel 
into (out of) the horizon (localized at — > —00) as the time increases. The solutions tp± 
can be expressed in terms of the normalizable and non-normalizable solutions: 



(2.20) 



where A(±) and are connection coefficients associated to Eq. ( 2.14 ), and determined 
as functions of w and q by imposing the boundary conditions on the field solution. 

The Schrodinger equation ( p. 14 ) can be expanded near the event horizon leading to 
the following expansion for ip^ and '0_|_: 



0± 



it iojr* 



1 + ai(±) ( 1 - — ) + a2(±) 



z 

1 

Zh 



+ 



(2.21) 
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where the elhpses denote higher powers of (1 — z/ z^) and the coefficients appearing in the 
foregoing expression are 



ai(±) 



6 + (g^ + 4c)4 
4 lb 2iujZh 



\—^{ [22 ± iiiOZh + (q^ + 4c)zl] ai(±) - 2zl{q^ + 8c) - 9 + 40^2^}. 



(2.22) 



16 lb Aiuozh 



The ingoing and outgoing wave functions form a basis for any other wave function. In 
particular, the normahzable and non-normaUzable solutions ipi and ip2 can be expressed as 

ip2 = C(2)ip- + I'(2)^+ , -ipi = C(i)V'- + ^5(1)^+ • (2.23) 
From (l2l2C| ) and (|]2|) we find a useful relation between the connection coefficients 



C(2) I'{2) 



(2.24) 



At zero temperature, it is not possible to impose the ingoing (outgoing) condition at 
the "horizon" z = oo due to the presence of the background dilaton field. Instead, we may 
impose the regularity condition 

limV'o = 0, (2.25) 



where V'o is a solution of Eq. ( 2.14 ) at T = 0. This wave function can be decomposed as 

V'o = Ao)V'2+^(o)V'i, (2.26) 



where ^(o) and B^o) are the zero-temperature connection coefficients associated to (2.14). 



2.4 An analysis of the effective potential 



We analyze here the effective potential V{z), given explicitely by ( ^.16 ), as a function of r^,. 
At zero temperature {zh oo), the tortoise coordinate r^, reduces to —z and the potential 
takes the form 



1 



^ + {q' + 2c)z' + c'z' 



(2.27) 



which is the potential considered in Ref. [19], and that leads to the glueball spectrum of 
Eq. (p.4|). In a similar way, the limit c — > reduces the effective potential to 



V^.=o (z) 



fiz) 



fz^ + — + 



15 9^4 



4 ■ 4z^ 



(2.28) 



which is the potential obtained in Refs. [31-33] in the absence of the dilaton. The effective 
potential V{z) suffers a very interesting transition in shape when going from high to low 
temperatures. Such a transition is discussed below. 
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Figure 2: The behavior of the effective potential V for high- and intermediate-temperature values. 
High temperatures 

In Fig. |2|, we show the potential as a function of the tortoise coordinate in the high- and 
intermediate-temperature regimes for q = 0. We observe the absence of any potential well 
so that there are no bound states in these regimes. For any fixed value of the tortoise coor- 
dinate r^, we note that the potencial increases with decreasing temperature. However, the 
shape of the potential does not change as one goes from high- to intermediate-temperature 
values. Therefore, on basis of the behavior of V, we do not expect a significative change in 
the spectrum of quasinormal modes as a function of temperature, for all > 1, in relation 
to the <I> = case. 
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Figure 3: The behavior of the potential for low temperature values and q = 0. The horizontal 
dotted lines correspond to the masses of the glueball states at zero temperature. 
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Low Temperatures 

At very low temperatures (T^ <^ 1), the effective potential begins to develop a well near 
the AdS boundary. This is illustrated in figure |3| where we plot the potential for some 
selected values of temperature and q = 0. As the temperature decreases the potential well 
gets deeper until it reaches the zero-temperature form. In the regime of low temperatures, 



the potential (2.16) with vanishing wavenumber can be approximated by 



----, + 2c.^+c^. 



.2 „4 



2c2^ + 0^ 



(2.29) 



where we have considered (z/z/i)^ to be small, keeping only first-order terms in this pa- 
rameter. The foregoing expression can be written in terms of the temperature as 



V 



15 
4^ 



2c c^z^ 



(2.30) 



where one can see explicitly the zero-temperature potential (2.27) plus low-temperature 
corrections, near the boundary. It is interesting to represent this potential in terms of the 
tortoise coordinate r^,. Near the boundary and for low temperatures, one can approximate 
the function 2;(r^,) by 



1 



^4 



Then, the potential ( 2.30|) reads 



V 



2c + cV^ + 



15 
4r? 



12 



2 4 6 T-i4 
-C TT r^T . 



(2.31) 



(2.32) 



The above potential is a combination of a constant term 2c, an oscilator-like contribution 
proportional to r^, an infinity barrier at = 0, given by the term 15/4r^, and temperature 
dependent corrections. The negative sign of the thermal corrections, in the low-temperature 
regime, imply that the real part of the QN frequency uj will decrease with the temperature. 
Indeed, we will observe this behaviour in the study of the quasinormal modes in section |^. 

An important consequence of the presence of the potential well at very low tempera- 
tures is the appearence of trapped modes. These modes are associated with scalar glueballs 
and will show up as peaks in the spectral function (the imaginary part of the Green's func- 
tion) as will be discussed in the next section. 



3. The retarded Green's function 



In order to calculate the glueball retarded Green's function we follow the prescription of 
Ref. [5]. First we use the equation of motion ( |2.10D to write down the on-shell version of 
the scalar field action (^^ ): 

Sonsheii = -n- / d^xy/^e'"^ g"""" <l)dz(p\,^o , (3.1) 
4/^10 J 
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where we have disregarded the boundary term at the horizon z = z^. The Fourier trans- 
formed on-sheh scalar field (^(z, k) is decomposed as 

4>{z,k)=ct>k{z)(l,o{k), (3.2) 

where (j)k{z) is the "bulk to boundary propagator" subjected to the boundary condition 

lim</)fc(z) = 1. (3.3) 

2:^0 

According to the Minkowkskian prescription of Ref. [5] , when computing the retarded 
Green's function one must consider a propagator (pki^) which satisfy an incoming- wave 
condition at horizon in the finite-temperature model, 



lim e ^ (-r 



Z \ -3/2 



(pkiz) 



1 



or a regularity condition^ at infinity {z — > oo) in the zero-temperature case, 



lim e 2 

z—>-oo 



/ Z\ -3/2 

Using the decomposition ( |3.2D , the action (|3.1| ) takes the form 

d'^k 



S. 



on shell 



(2vr)4 

The retarded Green's function is then given by 



-k)T{k,z) (po^k) l^^o . 



G^{k) = -2F{k,z) 



o hm ^ 



-gg^'e-''<l)l{z)dAk{z). 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Since 4>k{z) is a solution of Eq. ( |2.12D satisfying the normalization condition ( |3.3D and 
additionally the boundary condition (3.4) in the finite-temperature case and ( |3.5| ) in the 
zero-temperature case, it can expressed in terms of the wave functions defined in the last 
section as 

(3.8) 



, , N 3/2 £|i V'a 



{a) 



where the index a takes the value for the zero temperature case and — for finite temper- 
ature. Using the expansions (p.l7 ) and ( 2.18 ) for V'l and ^^2) we find 



/V2 



(c-h2a2i)e~^ -F461n(ce^) + 2 ib + a2ic + alA + ARe 



(a) 



A 



ilm 



(a) 



(a) 



A, 



(a) 



where e is an ultraviolet regulator. In the above equation we have used the relation 

24q ~ 87r2 ■ 



(3.9) 



(3.10) 



^The zero temperature condition (3.5) is equivalent to imposing the regularity condition lim^^oo 4>kiz) = 
0, for space-like momentum (fc^ > 0), and then performing a Wick rotation, similar to what was considered 
in Ref. [5]. 
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The first two terms in Eq. ( |3.9D are ultraviolet divergent and can be removed using a 
holographic renormalization procedure like that of Ref. [34]. It is interesting to observe that 
the imaginary part of the retarded Green's function at zero (finite) temperature depends 
only on the quocient i3(o)/^(o) (S(_)/^(_)). This quantity defines the well-known spectral 
function which has been calculated in Ref. [35] for vector mesons in the soft-wall model, 
and for other fields and holographic models, for example, in Refs. [26,36-39]. 

3.1 The zero-temperature case 

In the zero temperature case, the bulk to boundary propagator has an analytical solution 
consisting on a combination of Kummer and Tricomi confluent hypergeometric functions. 
The Kummer function is ruled out by the regularity condition leading to 



(3.11) 



where k'^ = —u'^+q^ and lA{a, 6, w) is the Tricomi function. Prom the asymptotic expansion 
of the Tricomi function we find near the boundary: 



2 Z 2 (j)^[z) 



_3 
Z 2 



|l -h a2iz^ + bz"^ ln(cz^ 



+ 



1 ^_^(^ ^ 
8^ 8c 8c Uc ^ 





1 




L4c J 



*(1) - ^-(3) 



2 4 
C Z 



+ 0{z^) + 0{zHncz^)y 



(3.12) 



where ^(x) is the digamma function and a2i and b are given in Eq. ( |2.19| ). From the above 
expansion we can extract 



(0) 



A, 



(0) 



+ 



8c ~ 8c V4^ 



+ l)[*(— + 2) -*(l)-^(3)]}c^ (3.13) 
— n, it is possible to extract an imagi- 



Since the digamma function ^{x) has poles at x = 
nary part of the retarded Green's function: 

ImG^ (A.) = -^Im^ = ^ ^ f ^ + l) limIm^f- + 2 

^=0^ ' 87r2 Am 87r2 8c V4c / e^o V4c 



HO) 

8^ 8c V4^ 



1 / 
^n-Fl V 

n=0 



4c 



+ 2 + n 



le 



(3.14) 



Hence, the spectral function at zero temperature is a sum of delta functions localized at 
/c^ = — 4c(n -|- 2). These poles are interpreted as the glueball masses at zero temperature. 
This result is consistent with the Feynman correlator obtained in Ref. [40]. 

3.2 Spectral function at finite temperature 
3.2.1 General procedure 

In order to find the spectral function at finite temperature we must calculate the imaginary 
part of the retarded Green's function. We see from Eq. ( |3.9| ) that ImG^(A;) depends on 
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the ratio a quantity that can be calculated by a numerical procedure consisting 

in the steps described below. 

Firstly, we solve numerically the Schrodinger equation ( p. 14 ) for a wave function that 
tends asymptotically to the normalizable solution ipi{z) as z ^ 0. In other words, we 
integrate numerically this equation from a point z = Zi near the AdS boundary to a point 
z = Zf near the horizon, using as "initial" conditions the truncated expansions 

= (1 + auzf) , d.M^^) = ^ zf/' (5 + 9anzf) . (3.15) 

The same numerical integration procedure is employed for a wave function that tends to 
the non-normalizable solution ip2iz) as z — > 0. In this case, we use the following conditions 
near the boundary: 

'>p2{zi) = z'^^"^ (l + a2izf) +b'4)i{zi) In (czf ), 

1 „5/2 2b (3.16) 

dzi'2{zi) = 7j + "■2iZi] +bdzipi{zi) In (czf) H i^i{zi). 

z Zi 

The expressions ( FTsI) and (III) are obtained from the near boundary asymptotic behavior 
for -01 and -02 discussed in the last section. 

With the numerical results for ■01 and '(/'2 at -2 = zj we construct the vectors 



(i = l,2), (3.17) 

with tp±{zf) and dzip±{zf) given by the near-horizon expansions ( 2.2l| ) truncated at second 
order in (1 — Zf/zh). Then we extract the coefficients Cy) and by inverting ( 3.17 ): 

(i = l,2). (3.18) 



T^U)) \dzi^-{zf) dz^+{zf)j \dzi^^,){zf)j 



Finally, from Eqs. (|]2|) and ( p^ we obtain 

■S(-) ^ ^ dz'4'-{zf)ilj2{zf) - i;-{zf)dz'il^2{zf) 

A-) ^(1) dzip-{zf)ilJi{zf) - ip-{zf)dz4^i{zf)' 



(3.19) 



This way we can obtain numerically using the numerical results for ipi, ip2, dzipi 

and dzip2 at a point z = Zf near the horizon. 

3.2.2 Numerical results 

The imaginary part of the retarded Green's function gives us the spectral function 

n{uj,q) = -21mG^{u;,q) = ^ Im-^ . (3.20) 

We present here some numerical results for the spectral function obtained using the proce- 
dure discussed above. We show in figures ^ and |5| the dependence of the spectral function 
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with the frequency uj = uj/\/c for selected values of temperature T = 7rT/-y/c and mo- 
mentum q = qj ^/c- We can see the presence of various peaks that correspond to the 
poles of the retarded Green's function. As commented above, these poles are associated 
with the frequencies of the black-hole quasinormal modes and are interpreted as glueball 
states in the dual gauge theory. Indeed, near one of the peaks the spectral function can be 
approximated by a Lorentzian function: 

n^.q)- ,~ ~\^~2 ' (3-21) 
[to -OJrY +UJf 

where Wjj and ujj are respectively the real part and the negative of the imaginary part of 
the black-hole quasinormal frequencies. The half- width of the distribution ( 3.21| ), given by 



uii, is interpreted in the dual gauge theory as the inverse of the glueball lifetime (multiplied 
by the constant \/c). The real part 5^, for g = 0, gives a measure of the glueball mass at 
finite temperature. 




Figure 4: Spectral functions for g = and selected values of temperature in units of N'^/A-k^. 
Note that the first peak of the spectral functions for the four cases almost coincide around = 8. 

We consider in Fig. ^ the zero-momentum case for various values of temperature. Note 
that as the temperature increases the number of peaks in the spectral function decreases 
while their width increase. Accordingly, the number of glueball quasiparticle states and 
their lifetimes decrease with the temperature. If we continue increasing the temperature 
we will reach a phase (for > 0.1) in which there are no more peaks. This suggests a kind 
of glueball melting in the sense that the glueball excitations disappear. A similar result was 
recently obtained in Ref. [35] when studying finite-temperature effects on the spectrum 
of vector mesons in the soft-wall model. It is important to remark that for the range of 
temperatures in which we found glueball excitations (peaks in the spectral function), the 
theory is in an unstable phase. 
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It is interesting to analyze the spectral function of Fig. ^ on basis of the behavior of the 
effective potential of Fig. ^. For instance, for = 0.02 there are three modes "trapped" in 
the potential well, which are represented by dotted lines. These quasi-stationary trapped 
modes correspond to the three sharp peaks that appear in the spectral function TZ{u!,0). 
Analogously, for = 0.03, the two sharp peaks at cD'^ ~ 8 and 5^ ~ 11 are associated to 
the modes trapped in the potential well. 

Note that the dotted lines shown in the effective potential graph of Fig. ^ correspond 
to the zero temperature modes. However, as shown in Fig. ^ the masses of the glueballs, 
given by the positions of the peaks, vary with the temperature. Nevertheless, this effect 
does not change the qualitative analysis made above. 

Another important feature shown in Fig. ^is that, for a fixed temperature, the width 
of the peaks increases with the frequency. This can be explained in terms of the tunneling 
effect of the scalar-field waves in the potential barriers of Fig. |3|, which is enhanced when 
the frequency becomes comparable with the height of the potential well. Then, the glueball 
states at finite temperature become less stable as the energy increases. 

14000 
12000 
10000 
O 8000 
6000 
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2000 


6 7 8 9 10 11 12 13 14 15 
co^/c 

Figure 5: Spectral function for = 0.05 and selected values of momentum q in units of N^/Att'^. 

In figure ^, we plot the spectral function for a fixed value of temperature = 0.05 
and various values of wavenumber q. In this picture one can see that the position and the 
width of the peak increases with the momentum. This result indicates that the glueball 
energy increases with the momentum, as expected, while the glueball lifetime decreases. 

In this section the frequency oj was taken as a real quantity for the computation of the 
spectral function TZ{u!,q). However, when studying the black-hole quasinormal modes (the 
poles of G^) in the next section, we will look for frequencies of the form ujqjs, = ujjt — iujj 
with (jjj and being (positive) real quantities. 
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4. Quasinormal modes and the scalar glueballs at finite T 



Quasinormal modes (QNMs) are solutions for classical field perturbations on a black hole 
geometry (see Refs. [9,41,42] for general reviews on the subject). These solutions are 
subjected to an absorption condition at the horizon and a regularity condition at the 
spacetime boundary in the AdS black hole case. The absorption condition leads to a 
purely ingoing wave at the horizon while the regularity at the boundary translates into a 
Dirichlet condition for tp, for the case of a massless scalar field discussed here. 

The analysis of the effective potential in section || gives additional support to the above 
quoted boundary conditions. The fact that the effective potential diverges at the spacetime 
boundary z = makes it natural to impose a Dirichlet condition, while the vanishing of 
the potential at z = Zh and the one-way membrane property of the horizon lead naturally 
to ingoing plane wave solutions. 

The absorption condition can be implemented by means of the solution ^_ (z) of section 
^ which behaves like an ingoing wave at the horizon, according to Eq. ( 2.211) . Imposing 
Dirichlet boundary condition on ip-{z) at z = implies that the connection coefficient 
A(-){uj, q) vanishes. Since the retarded Green's function ( ^ ) diverges when A(-){uj, q) = 0, 
we conclude that the black-hole quasinormal modes correspond to the (complex) poles of 
the glueball retarded Green's function. 

In the following, we will show the numerical results obtained for the complex roots 
ljqn = ojji — iuJi of the equation A(^) {oj,q) = 0. These complex roots are the frequencies of 
the scalar-field quasinormal modes of the black hole in the presence of the soft- wall dilaton 
background. Previous studies of quasinormal modes in the context of gauge/gravity duality 
can be found, for example, in Refs. [27,43-55]. 



4.1 Methods 

In the last decades, different methods have been developed to numerically calculate the 
quasinormal modes of black holes. We employ here two of these techniques to find out the 
QNMs of a scalar field in the soft-wall model. The first one is the traditional Frobenius 
power series method, following the approach of Horowitz and Hubeny [27]. In this case 
the problem of finding quasinormal frequencies is reduced to that of finding the roots 
of a polynomial equation. The second method is based on ideas from the study of the 
ressonant scattering problem in atomic and nuclear physics. It has been recently applied 
to the calculation of QNMs [28] and is presently known as the Breit-Wigner resonance 
method. 



4.1.1 The power series method 

The power series method of Ref. [27] is of particular interest for the present work since 
it is suited for asymptotically AdS spacetimes. We start with the scalar field equation 
of motion written as a Schrodinger-like equation ( ^.14 ). As discussed above, QNMs are 
defined as solutions which are purely ingoing waves at horizon, ifj — > e"*'^''*, and satisfy a 
Dirichlet condition at z = 0. In the power series method, it is convenient to introduce a 
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new function = e^^^^*4>- Then, Eq. ( 2.14| ) becomes 



The above equation has a regular singular point {z = z^) in the region of interest, < z < 



Zh- So, according to the Fuchs's theorem, at least one of the solutions of Eq. can be 



written as a Frobenius series [56]. In fact, from the solution ( 2.21| ) for ij)- one can see that 



the corresponding function can be written as a power series of the form: 

oo / \m 

(^_(z) = J^a„(_) 1-- . (4.2) 

m=0 ^ 



Substituting Eq. (|4.2| ) into (4J.) we obtain a recurrence relation for the coefficients a^'s. 
Since Eq. ([4.lD is linear we can set oq = 1 and calculate this way any coefficient am(-) as 
a function of the frequency oj = w/y^, the wavenumber q = q/\/c, and the temperature 
T = ttT / y/c. The coefficients ai(_) and a2(_) are given explicitly in equation ( p. 22 ). Then, 
imposing the Dirichlet condition at AdS boundary, we obtain 

oo 

(/p_(0) = ^a^(_, =0. (4.3) 

m=0 

Thus, the problem of finding quasinormal (QN) frequencies has been reduced to that of 
obtaining the roots of the equation ( [4.31 ). Since we cannot carry out numerically the full 



sum in expression (4.3), we truncate the series after a large number of terms and look for 
the zeros of the resulting partial sum. The precision of the results is then verified through 
the relative variation between the roots of two successive partial sums. For each pair of q 
and T, we find a set of QN frequencies tD, labeled by the principal quantum number n and 
ordered in a set beginning by the roots with the lowest absolute values for the imaginary 
part. 

4.1.2 The Breit-Wigner method 

Although the power series method is suitable to find QN frequencies in the intermediate- 
and high-temperature regimes (T > 1), its convergence gets weak at very low temperatures 
(T <C 1). Fortunately, in this regime the behavior of the effective potential V{r^,) enable 
us to use an alternative method based on the Breit-Wigner theory of resonances. 

The theory of Breit-Wigner resonances was originally employed in gravitational prob- 
lems to find the resonant frequencies associated with ultra-compact relativistic stars [57- 
59], and only recently this technique was used to find QN frequencies of black holes, par- 
ticularly of very small Schwarzschild-AdS black holes [9,28]. An important element for 
the efficiency of this method is the form of the effective potential: V{r^:) should have the 
form of a very deep well, so that quasi-stationary "trapped" states can exist and the scalar 
field waves can only leak out to the horizon by tunneling through the potential barrier. 
This is exactly the situation found in section |^ when analyzing the effective potential for 
a massless scalar field in the soft-wall black hole geometry. Below we present the general 
idea of this method in the form used in asymptoticaly AdS spacetimes. 
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The resonance method deals with normahzable solutions, which corresponds here to 



choosing the solution ipi of Eq. (2.17). As discussed in section g, the function ipi can be 
written as a combination of ingoing and outgoing waves tlj±, as shown in Eq. ( |2.23| ). This 
equation determines the asymptotic behavior of tpi at the horizon: 

tj;^ = Cie"^'^''* + Vie''^''* = acoswr* + /Jsinwr,, (4.4) 

where a = Ci +Pi and (5 = —i{Ci —Di). The QNMs are normahzable solutions with com- 
plex frequencies loq^ = uJr — iui such that there are no outgoing waves. This corresponds 
to set Pi = 0. Analogously, the complex conjugate of a QNM corresponds to the vanishing 
of the ingoing coefficient (Ci = 0). Assuming that the QNMs are simple zeros of Pi we 
expect that near these modes Pi ~ (w — lor) + iuji and Ci ~ (a; — lOh) — iui. Then it turns 
convenient to introduce the real quantity 

+ /?2 = iCiVi const, x [{lu - LOnf + u]] . (4.5) 

The Breit-Wigner resonance method is based on the minimization of this quantity, consid- 
ering a; as a real variable [57,59]. Note that this is far simpler than solving the complex 
equation Pi = when using the numerical integration procedure of section ^. 

The quantity + will present sharp dips at lo = corresponding to the real parts 
of the QN frequencies. Once we have located a minimum of -|- on the line of real 
uj at some uj = ujn, the imaginary part of the QN frequency iVj is obtained by means of 



a parabolic fit of the Breit-Wigner formula (4^). To set the accuracy of our numerical 



results, we use the alternative (and equivalent) expressions for lj^: 

= -p/a' = a/(i', (4.6) 

where a prime indicates a derivative with respect to w^, evaluated at the minimum [28,58]. 
The Breit-Wigner method has, in general, a good convergence as long as the condition 
ujj <^ LOn is satisfied. 

It is interesting to notice the relation between the dips of -|- and the peaks of the 
spectral function TZ{uj, q), discussed in section |^. Each minimum of (|4.5| ) appears as a sharp 
maximum of the spectral function. This connection can be understood by considering the 
approximate expression 

7^ = -2ImG^ ~ (A-))"^ ~ (^i)"^ ■ (4-7) 



Therefore, the simple pole hypothesis for the retarded Green's function that leads to ( 3.21| ) 
is directly related to the simple zero hypothesis for the coefficient Pi. 

4.2 Numerical results for the QN frequencies 
4.2.1 Temperature dependence 

Here we present our numerical results for the QN frequencies when varying the temperature 
while keeping the momentum q = 0. The power series method can be applied to a wide 
range of values from high temperatures to ~ 0.022. At lower temperatures, the power 
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series method has convergence problems but, as we discussed above, it is precisely in this 
regime that the wells of the potential are deep and the dips in c? + are sharp enough 
to make the Breit-Wigner method efficient. We have used the Breit-Wigner method for 
< 0.05. A comparison of the results obtained with the two methods is presented in 
Table [l|. The approaches are complementary, and show a very good agreement for the 
range 0.025 < f ^ < 0.04 . 
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Breit-Wigner 








^? 
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^? 




0.0493827 


7.52779 


1.41127 X 10" 


-4 


7.53270 


1.43646 X 10" 


-4 


0.04 


7.72086 


1.31516 X 10" 


-6 


7.72092 


1.31549 X 10" 


-6 


0.0330579 


7.82420 


2.12446 X 10" 


-9 


7.82420 


2.12446 X 10" 


-9 


0.0277778 


7.88073 


8.61149 X 10" 


13 


7.88073 


8.61151 X 10" 


13 


0.0236686 


7.91529 


1.04760 X 10" 


16 


7.91529 


1.04873 X 10" 


16 


0.0204082 








7.93787 


2.23960 X 10" 


20 



Table 1: A comparison between the values obtained for the fundamental [n = 0) QN frequency 
computed using the power series and the Breit-Wigner resonance methods. 

The values of Wj shown in Table ^ for the Breit-Wigner method were calculated using 

3) for -|- The alternative expressions (4.6) were used to 



the parabolic fit of Eq. 
check the consistency of the results in the range of temperature values where the power 
series method does not converge. In this case, the number of significant figures retained in 
uJi was chosen so that all the alternative ways of determining ujj give the same value. 





Figure 6: Numerical results for the square of the real and imaginary parts of the QN frequencies, 
Lu'jf/c and w^/c, for the first five quasinormal modes n = 0, 1, 4, with q = 0. 

In Fig. ^ we show our results for the square of the real and imaginary parts of the 
QN frequencies as functions of T'^ for the first five modes n = 0, 1, 4 with g = 0. We 
observe that, in the zero temperature limit, the real part of the QN frequencies coincide 
with the corresponding glueball mass spectrum, given by Eq. (p.4|). Then, increasing 
the temperature tD^ decrease until they reach a minimum value (tD^)™'*" at some critical 
temperature, different for each mode. The value of the critical temperature is lower for 
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higher overtone numbers n. For higher temperatures, cj^ increases approaching a Unear 
dependence on the temperature. On the other hand, we observe in Fig. ^ that (the negative 
of) the imaginary part of the frequencies lji increase monotonicaUy with the temperature 
approaching constant slopes. In Table ^, we show the values of the minimum of the curve 
tD^ versus T-^ for the first five quasinormal modes. We also show there the numerical values 
for uj] and the quality factor Q = lor/2loj evaluated at the critical temperatures. Notice 
that the values of Q in the minimum of ujn are such that Q ^ 4 and so, for temperatures 
equal or higher than the critical value, there are no more quasiparticle peaks corresponding 
to glueballs in the spectral function 71{lo,0). Quasiparticle states show up for higher values 
of Q. 
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/~2 \min 




Q 





0.122069 


6.85489 


0.108102 


3.98156 


1 


0.0637774 


9.92861 


0.166347 


3.86284 


2 


0.0444836 


13.0406 


0.241963 


3.67066 


3 


0.0343047 


16.1605 


0.319899 


3.55378 


4 


0.0279468 


19.2844 


0.398241 


3.47936 



Table 2: The minimum value of the real part of QN frequency for the first five modes, n = 0, 1, 4, 
and the corresponding values of cD^, T^, and the quality factor Q. 



The decreasing of the real part of the quasinormal frequencies at very low temper- 
atures, shown in figure 0, is related to the fact that the very low temperature potential 



given by ( 2.32 ) behaves as a harmonic oscillator with an infinite barrier r^^ and a tem- 
perature dependent negative correction —T^r^. This negative correction is responsible for 
the observed decrease in u)\ for < 0.1 . 

The exact way the functions a;^(r^) and u}f{T'^) deviate from the linear behavior, 
observed in high temperatures, varies with the level n. Within a certain accuracy, the fre- 
quency of each mode deviates significantly from a straight line only in the low temperature 
regime. For instance, this deviation for the fundamental mode occurs below = 0.25 . For 
higher overtones, as n increases, the deviations from the linear behavior occur at decreasing 
temperatures. 

We now compare our results for the soft-wall model with the pure AdS black-hole case, 
which corresponds to the limit $ = 0. Some results for the fundamental mode (n = 0) 
obtained with the Horowitz-Hubeny power series method are shown in Table ^ for q = Q 
and selected values of temperature T. We introduce for the real and imaginary parts of 
the QN frequencies the fractional differences 

^ fl^ ^ RJ± ^ 

where Vo = ljo/2'kT is the normalized frequency commonly used in the literature [36,47,48] 
and the superscript $ = indicates the results for the AdS black hole without dilaton. 
Since the frequencies tt)*"° and trf"" are constants independent of the temperature, for a 
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fixed level n the fractional differences A^j varies with T due to the changes in the values 
of =U}nj/2f. 

As expected, one can see a significant difference in the numerical results for the soft- 
wall and pure AdS cases for low temperatures. For higher values of T, the effect of the 
background dilaton field is small, and the QNMs are essentially the same as those of the pure 
black hole solution [46,47]. This is related to the fact that the effective potential V{r^) is a 
monotonic function without wells both in the intermediate- and high-temperature regimes. 



f2 


ul 




Ah 


A, 


Q 


100 


976.020 


750.874 


0.00150017 


0.0023548 


0.570054 


25 


246.247 


185.138 


0.0060918 


0.00923686 


0.576645 


4 


42.1859 


27.1783 


0.0410586 


0.0509827 


0.622934 


1 


13.6137 


5.45053 


0.182798 


0.150013 


0.790204 


0.25 


7.40192 


0.672991 


0.744311 


0.402652 


1.65820 


0.04 


7.72086 


1.31516 X 10"^ 


3.45374 


0.997912 


1211.47 



Table 3: The frequencies of the fundamental QNM for some selected values of temperature T, 
calculated with q = 0. We also show the fractional differences in relation to the frequency values 
in the absence of dilaton, and Aj^, and the quality factor of the oscillations, Q. 



4.2.2 Dispersion relations 

In this section we investigate the momentum dependence of the scalar-field black hole 
quasinormal modes in the soft-wall model. In Table we list some of our numerical results 
for the real and (negative of) imaginary parts of the QN frequencies, ujR{q) and uji(q), for 
selected values of T and two different momenta, q = 1, 10. We also compare uj{i{q) with 
the standard relativistic dispersion relation to = \J (wJJ)^ + q^^ where = ^^(g)!^^^ can be 
regarded as the glueball mass at finite temperature. We introduce the fractional difference 



KM_^M, ,4.3) 

which is a measure of the fractional difference between the numerical values obtained for 
ijjji and the relativistic energy-momentum relation for a massive particle. 

In the ranges of momentum and temperature considered here we have obtained very 
small values for Ag. The maximum value found for this quantity corresponds approximately 
to 1.3%, and appears for g = 10 and T-^ = 0.2. The overall behavior of the dispersion 
relation ^^(g) is such that A^ increases with the momentum and presents an oscillatory 
dependence in T for low temperatures, while increases monotonically for higher T . 

We present in Fig. ^ the variation of the real and (negative of) imaginary parts of the 
QN frequency with the momentum for selected values of the temperature. In the left panel 
of this figure, we illustrate the behaviour of the real part of the frequency in the range 
< g < 1. From this picture, one can see that 0Jii{q) increases with the momentum. In 
our study we found this behaviour for all temperatures investigated. 
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A. 




UJl 


A, 


0.2 


2.85042 


0.644367 


0.0010234 


10.4808 


0.562306 


0.0127022 


0.15 


2.80845 


0.452717 


0.000665236 


10.4244 


0.449494 


0.00824847 


0.1 


2.80271 


0.231549 


0.0024438 


10.3660 


0.317525 


0.00262035 


0.05 


2.91009 


0.0157488 


0.00267834 


10.3162 


0.143700 


0.00508183 



Table 4: Some data for the fundamental QN frequency of the gluebaU fluctuations for q — 1, 10 
and selected values of temperature. 
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Figure 7: The dispersion relations of the fundamental quasinormal mode n = for some selected 
values of temperature. 
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0.20 





2.66613 


0.647625 


1.54410 


0.15 


5.1149 


5.77971 


0.475506 


2.10302 


0.10 


8.5667 


8.9751 


0.319018 


3.13462 



Table 5: The maximum absolute value of the imaginary part of the fundamental QN frequency as 
a function of q, the corresponding values of Wh, T^, and the characteristic damping time r — l/uj, 
in miits of c^^/^. The numerical methods used here present bad convergence for < 0.05 . 

In the right panel of Fig. 0, we show the dependence of (the negative of) the imaginary 
part of the frequencies with the momentum. The power series method converged well 
for momenta in the range < g < 40 and = 0.2, but only for < g < 20 when 
= 0.05 . We found that the absolute value of the imaginary part of the frequency presents 
a maximum value, which means a minimum for r = l/wj, for some finite momentum. For 
temperatures above a critical value ~ 0.16, the maximum value of uji is achieved at 
g = 0. For temperatures below this critical value, uij increases with the momentum for very 
low g 's, reaches a maximum value cD™"^ and then decreases monotonically. Our numerical 
results suggest that lOj vanishes for infinitely large q, a behavior similar to that of the pure 
AdS black hole case investigated in Refs. [60,61]. 

In Table |5|, for some selected temperatures T, we present the maximum values of u)j and 
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the corresponding values of g, cj^j, and the thermahzation timescalc r = l/cD/""^. Below 
the temperature = 0.05 the convergence of our numerical methods gets weak, and we 
could not determine the value of w/""^ for < 0.05 . In the AdS black hole space without 
the dilaton, the (negative of the) imaginary part of the frequency is a decreasing function 
of the momentum q [46,47], which is in accordance with the behavior of uii for > 0.16 . 

5. Final comments and conclusion 

In this article we investigated the spectrum of scalar glueballs at finite temperature in the 
soft-wall model. Such a study was performed considering a scalar field in an AdS black 
hole spacetime in the presence of a dilaton field as the gravity dual of finite temperature 
scalar glueballs. This background corresponds to the deconfined plasma phase which is 
thermodynamically favoured at high temperatures. We considered here the AdS black 
hole space for all temperatures. For intermediate temperatures it represents a supercooled 
metastable (positive specific heat) plasma and for low temperatures an unstable (negative 
specific heat) plasma. 

We obtained the thermal and momentum dependences of the frequencies of the black- 
hole quasinormal modes. For g = 0, the real part of the QN frequencies furnishes the 
glueball masses for low temperatures. In the zero temperature limit our results give the 
glueball spectrum found in Ref. [19] for the soft-wall model. Surprisingly, for very low 
temperatures and low fixed momentum, these masses decrease with increasing temperature. 
After reaching a minimum value, the real part of the QN frequencies grows monotonically 
with the temperature, as in the case without the dilaton field. 

At zero temperature there is an infinite number of stable glueball states corresponding 
to delta function peaks in the spectral function. As the temperature increases, however, 
the peaks spread, and for a fixed low temperature and fixed low momentum, there is only 
a finite number of sharp peaks, since the width of the peaks increase with the frequency. 
The sharp peaks correspond to QNMs whose complex frequencies have small imaginary 
parts. These modes can be interpreted as glueball quasiparticle states with decay rates 
inversely proportional to the imaginary part of the frequencies. This imaginary part is also 
proportional to the width of the peak of the spectral function. 

We studied the dispersion relations for the glueballs in the soft wall at finite temper- 
ature. This way we found the effect of the quasiparticle motion on its lifetime inside the 
plasma. The connection of the complex QN frequencies with the mass and the decaying 
time of quasiparticles in the dual plasma gives an important tool to explore both sides of 
the duality using known elements from the gauge theory and the physics of black holes. 

The process of spreading of a peak as the temperature increases is interpreted as the 
melting of the quasiparticles in the thermal bath. We found that the scalar glueball states 
melt at temperatures much lower than the confinement/deconfinement temperature. In 
the soft-wall model the phase transition occurs at T = Tc with (7rTc)^/c pa 2.39, while we 
found that, at (7rT)^/c ~ 0.1 all the quasiparticles have melted. 

Our results indicate that the glueball formation in the black hole phase occurs at very 
low temperatures compared to the deconfinement temperature. So, from the point of view 
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of the soft-wall model, it would be difficult to observe experimentally the contribution of 
the black hole phase to glueball formation unless the plasma cools down in a very small time 
scale, compared to that of the thermal phase transition. Recall that in the soft-wall model 
there are two different thermal phases, corresponding to different backgrounds. As shown 
in ref. [22] there is a first order phase transition between them, implying that these two 
phases can coexist. For T < Tc, the thermal AdS space without a black hole is dominant. 
In this confining background the glueballs are formed at temperatures just below Tc having 
constant masses and zero widths. 

As already mentioned, the soft- wall model is not a solution of the 5D Einstein equa- 
tions since it does not take into account the dilaton backreaction. There are some recent 
AdS/QCD models at zero temperature that overcome this problem [62-64]. At finite tem- 
perature, the thermodynamics of AdS/QCD models with dilaton backreaction was consid- 
ered in [65,66]. It would be interesting to extend the analysis of the present work to these 
kind of models. This way one could distinguish the aspects that represent general char- 
acteristics of quasiparticle excitations in the plasma from those specific to the holographic 
soft-wall model considered here. 

After completing this work we noticed the very recent article [67] where glueballs and 
mesons in the soft-wall model at finite temperature are studied, with some overlap with 
this work. 
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